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基于 上 下 行 对 偶 的 CRN 下 行 功率 分 配 和 波束 赋 形 算法 
季 中 恒 ， 季 新 生 ， 黄 开 枝 


(国家 数字 交换 系统 工程 技术 研究 中 心 ， 郑 州 450002) 


摘 要 : 认 知 无 线 网 络 (CRN) 在 underlay 工作 模式 下 的 多 用 户 下 行 功率 分 配 和 波束 赋 形 问题 研究 中 存在 通用 的 SDR 
算法 计算 复杂 度 高 ， 实 用 性 受 限 以 及 优化 问题 中 忽视 主 网 络 (PN) 对 认 知 用 户 〈SU) 的 干扰 等 问题 。 针 对 这 些 问题 ， 
首先 在 CRN 网 络 模型 中 增添 PN 对 SU 的 干扰 ， 生 成 优化 问题 ; 而 后 基于 上 行 和 下 行 的 对 偶 特 性 ， 采 用 虚拟 功率 ， 将 
优化 问题 进行 形式 变换 ， 成 为 上 行 功率 分 配 和 波束 赋 形 问题 ; 得 到 能 够 简便 、 快 速 求解 的 迭代 算法 。 分 析 了 算法 的 收 
敛 特性 ， 得 到 了 收敛 条 件 ; 并 进一步 计算 了 算法 的 复杂 度 ， 结 果 表 明 优 于 SDR 算法 。 数 值 仿真 显示 ， 算 法 收敛 很 快 ; 
而 且 表 明 主 网 络 基站 (PBS) 发 送 功率 的 变化 影响 可 行 解 区 域 ; PBS 发 送 功率 的 增加 会 导致 CRN 下 行 功率 增 大 ， 影 响 
较 显 著 。 
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CRN downlink power allocation and beamforming algorithm 
based on duality of uplink-downlink 


JiZhongheng, Ji Xinsheng, Huang Kaizhi 
(National Digital Switching Systems Engineering & Technological Research Center, Zhengzhou 450002, China) 


Abstract: There are some problems on multi-user downlink power allocation and beamforming in a cognitive radio network 
(CRN) which works in the underlay mode. They are the high complexity of the conventional SDR algorithm, the limited 
applicability and ignoring the interferences of the primary network (PN) to the secondary user (SU) , etc. Aiming at these 
problems, this paper firstly added the term of interference of the PN to the SU to the CRN model and formulated the optimization 
problem. Then based on the duality of uplink-downlink, it obtained the simple and fast iterative algorithm by introducing virtual 
power, transforming the optimization problem into the problem of uplink power allocation and beamforming. And it analyzed 
the convergence of the algorithm to obtain the convergent conditions. It also reckoned the algorithm complexity. The result 
demonstrates that iterative algorithm is superior to the SDR algorithm. Numerical simulation results show that the algorithm 
converges faster, and that the variation of transmitting power of the primary base station (PBS) could affect the feasibility region. 
They also show that the increasing of transmitting power of the PBS could lead to the downlink power increasing. The influences 
are notable. 
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本 underlay。 在 interweave 工作 模式 , CRN 只 能 在 PN 不 工作 的 频 

9 点 上 接 入 ， 而 一 旦 PN 要 占用 该 频 点 ，CRN 必须 退出 ， 这 样 就 
认 知 无 线 网 络 (cognitive radio network,CRN) 技术 能 有 效 有 可 能 会 使 CRN 进行 多 次 的 接 入 和 退出 , 使 其 无 法 正常 工作 。 
地 利用 无 线 频 谱 ， 认 知 无 线 网 络 用 户 ， 即 次 级 用 户 (secondary ”在 overlay 工作 模式 ，CRN 需要 获知 PU 的 码 本 和 消息 知识 ， 

user，SU) 可 以 在 不 影响 授权 用 户 ， 即 主 用 户 (primary user， 通过 编码 方法 来 消除 对 PN 的 干扰; 但 是 这 一 条 件 ， 在 多 数 场 

PU) 的 前 提 下 , 接 入 频谱 。 这样 就 形成 了 两 个 网 络 , 一 个 是 PU  ” 合 下 ,并 不 能 得 到 满足 , 因而 CRN 会 影响 PN 的 工作 。 在 underlay 


组 成 的 主 网 络 (primary network，PN)， 男 一 个 是 SU 组 成 的 工作 模式 ，CRN 通过 信道 条 件 来 控制 自身 的 发 送 参 数 ， 实 现 与 
CRN， 两 个 网 络 在 无 线 频谱 空间 是 县 加 在 一 起 的 。 PN 共享 频谱 。 
CRN 有 三 种 工作 模式 ， 分 别 是 interweave，overlay 和 工作 于 underlay 模式 下 的 CRN 多 用 户 功率 分 配 和 波束 赋 
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当前 CRN 研究 的 一 个 重点 。 当 SU 和 PU 的 数目 较 少 
时 ， 可 采 / (zero forcing，ZF) 算法 ， 来 实现 CRN 和 PN 
的 频谱 共享 由 3。 但 是 当 SU 和 PU 的 数目 增多 时 ,ZF 算法 就 无 
需 通过 解 有 约束 条 件 的 最 优化 问题 来 实现 CRN 和 PN 
的 频谱 共享 64， 优 化 目标 是 下 行 功率 即 CRN 基站 (CBS) 的 
发 送 功率 最 小 , 约束 条 件 有 以 信 干 品 比 (SINR) 来 度量 的 CRN 
的 服务 质量 QoS (quality-of-service )、 波 束 赋 形 矩阵 及 CBS 对 


S 


PU 的 总 干扰 等 。 研究 中 存在 以 下 的 问题 : 一 是 最 优化 问题 求解 
的 通用 实现 算法 SDR (semi-definite relaxation) 有 较 高 的 计算 
复杂 度 钻 , 实用 性 受 限 ; 二 是 在 优化 问题 中 忽略 由 PN 产生 的 对 
SU 的 干扰 [5 5， 这 是 缘 于 PN 的 服务 优先 级 高 。 但 是 当 在 满足 
PN 服务 要 求 的 前 提 下 ， 如 果 不 把 PN 所 产生 的 对 SU 的 干扰 计 


入 CRN 的 工作 约束 条 件 中 , 将 使 CRN 的 QoS 条 件 不 严格 ， 导 
致 优化 结果 在 实际 使 用 时 产生 较 大 的 不 确定 性 。 三 是 虽然 提出 
了 非 SDP 算法 ,如 文献 [7] 的 迭代 算法 , 但 是 没有 给 出 明确 的 收 
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图 1 CRN 网 络 模型 


SU 的 接收 信号 y (7) : 


> A Us, (四 十 


=1,iz 


yD = Vpe he Ur sD+ 2 


JP h,. ss(D + zi(D) (1) 


敛 性 条 件 和 复杂 性 分 析 。 SO (j=12,..sK ) CBS 发 送 给 SU j 的 信号 ， 
对 以 上 的 问题 ， 对 工作 于 underl JJ CRN 2 SR 2 
本 有 人 lwol )=1 s,(1) 是 PBS 的 发 送信 号 ， Ellsol j=! Zl) 
行 功率 分 配 和 波束 赋 形 问题 作 进一步 的 研究 , 参照 文献 [11,12]， 
建立 CRN 网 络 模型 ， 在 约束 条 件 中 增添 PN 对 SU 的 干扰 项 ， 是 加 性 高 斯 白 噪声 ， z (1) ~ N(0,a2) 。 
得 到 优化 问题 ， 在 文责 I 基础 上 ， 基 行 对 偶 原 理 ， K 
得 到 优化 问题 ， 在 文献 [7,10] 的 基础 上 ， 基 于 上 下 行 对 侦 原 于 CBS 对 PU1(1=12..,L ) 的 干扰 信号 为 部 DMK ua 。 
引入 虚拟 功率 ， 将 下 行 功 率 分 配 和 波束 赋 形 问题 转换 为 上 行 虚 气 
拟 功率 分 配 和 波束 赋 形 问题 ， 先 得 到 上 行 虚拟 功率 和 波束 赋 形 SU 的 SINR 为 
向 量 ， 而 后 再 由 对 应 的 变换 得 到 下 行 功 率 ， 建 立 迭 代 关 系 ， 用 peul 
迭代 算法 快速 、 简 便 地 求解 优化 问题 。 并 对 算法 的 收敛 特性 和 SINR = < 万 内 2 一 
大 性 AD 六 全 一 人 日 A 于 九 性 话 旦 . > pilh ul +p,lhl + 

复杂 性 分 别 进行 了 分 析 ， 得 到 了 收敛 条 件 和 复杂 性 度量 值 。 k 

通过 数值 仿真 , 检验 了 迭代 算法 的 收敛 特性 ; 分 析 了 PN 基 pus hh u, 
站 (PBS) 的 发 送 功率 与 可 行 解 区 域 之 间 的 关系 ; 下 行 功率 与 pur hihe + pohihs + o? (2) 
PBS 的 发 送 功率 之 间 的 关系 等 。 仿 真 结果 显示 ， 迭 代 算 法 收敛 lk 
很 性 > 通过 过 几 才 I 次 失 失 代 后 即 可 收敛 ; PBS 发 送 功率 的 变 化 会 影 对 为 h, 的 维 数 为 I FP nl = hih, = hh 川 是 
响 可 行 解 区域 ， 功 率 增加 ， 可 行 解 区 域 减 小 PBS 的 发 送 功率 
对 下 行 功率 有 较 明 显 的 影响 ， 尤 其 是 变 大 时 ， 下 行 功率 的 变化 ”Euclidean 范 数 。 
更 加 显著 。 在 实际 应 用 中 ， 很 难得 到 信道 的 瞬时 特性 ， 一 般 是 使 用 期 

i 望 值 。 因 此 可 设 : SUk 与 CBS 之 间 的 信道 相关 和 矩阵 有 = 

1 CRN 网 络 模型 及 优化 问题 E(h.h*)，SUX 与 PBS 之 间 的 信道 相关 和 矩阵 R,,= E(h,.ha)， 
1.1 网 络 模型 用 RR 和 记分 别 蔡 换 hh” 和 hhs 可 得 : 

网 络 模型 如 图 1 所 示 ， 由 1 个 CBS 和 多 个 SU 组 成 CRN,， pus Ru, 
1 个 PBS 和 多 个 PU 组 成 PN。CBS 和 PBS 共享 频谱 。 实 际 网 SINR, = G3) 


络 中 PBS 的 数目 不 止 1 个 , 但 其 他 PBS 对 SU 的 干扰 远 没有 与 
CBS 共享 频谱 的 PBS 的 干扰 效应 显著 , 因此 对 分 析 结 果 几 乎 无 
影响 ， 可 以 略 去 。PU 的 数目 为 上 ，SU 的 数目 为 及，CBS 的 发 
送 天 线 数目 为 W ，PBS 为 单 天 线 发 送 ，PU 和 SU 均 是 单 天 线 
接收 。 

CBS 与 SU 之 间 的 信道 向 量 h = (hi,h 
的 功率 为 p=(p,p,,.…, pr) ， 波 束 赋 形 矩阵 y= (ww 
其 中 : Ur = (Ur Ui 20s Uw)” 。 

PBS 与 SU 之 间 的 信道 向 量 


,Ux) . 


H 
2， hn, ) ? 


h,。( 维 数 为 1), 发 送 的 功率 


K 
H 2 
>» Pi RU; + p,P, + ox 


i=L,i#k 


CBS 对 PU 1 (1=1,2,.…,L ) 的 干扰 信和 号 号 DA hi Us (CD ， 
i=1 


> hw | 本 -> Ur hehesd (4) 


可 设 :PU1 与 CBS 之 间 的 信道 相关 和 矩阵 PR ,= E(h hs)， 


替换 hh. ,hs 可 得 : 


Rr,, 天 +1 K+l 
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K 
L 二 > pu KilUi 。 
i=1 
关于 CRN 模型 的 两 点 说 明 : 
1) RR 有 以 下 的 特性 : 
(1) 忆 是 Hermitian 矩阵 ， 大 =12, 天 十 工 


由 忆 =E(h.h”) 可 得 : R: =E(hh’)”=E(hh!’)=R,。 


>0。 
2) 与 文献 [9] 将 PN 干扰 用 
型 引入 PBS 和 干扰， 能够 更 加 准确 地 


(5) 


(2) 忆 >0， 即 忆 是 正定 / 半 正 定 的 ，k =1,2,...,K+L。 


因为 : x*RXx = Xx!E(h.h’)x, = E(x hh’ x,) = 人 la 


随机 数值 表示 不 同 ，CRN 模 
省 述 SINR, ， 并 且 PBS 与 


SU 之 间 的 信道 的 时 变 特性 使 干扰 项 P PR， 对 SINR, 的 影响 显 
著 ; 进一步 影响 p、u, 和 等, 使 优化 问题 的 结果 更 符合 实际 。 


1.2 优化 问题 
下 行 功率 分 配 和 波束 赋 形 的 优化 问题 如 下 : 


K 
Pl1: Cs .7 
peur Ru 
St. SINR.= < .14; k=1,2,...,K 
* > piur RU 十 局 十 CO A 
i=1,i#k 
pr >=0;|u|=1; p, >0 
化 
1 = 2 pu Rd, < ;1=1,2,..,L 
i=1 K+l 
1 、 
XX 和 分 别 是 SINR, 和 了 的 门限 
K+l 
问题 P1 的 形式 转换 成 为 P2。 
K 
P2: on 7 
及 此 及 由 
St. SINR = < yr; k=1,2,...,K 
人 > Pi Ru;+ popP + Ox 
i=1,i#k 
pr >=0;|u|=1; p, =>0 


引入 松弛 变量 p<1 (1=1,2,.…,L), 作为 虚拟 的 PU 下 行 


功率 ， 可 得 问题 P3。 


K 
P3: Os 7 
piu¢ Ru, 
Sst. SINR,= 蕊 之 V3 天 三 1 2 天 
>》， pwr Ru, + p,P, + Ox 4 
i zk 
Pra 
SINR, ,= => Tal E 
0 2 pu Re, 6 


(6) 


(7) 


(8) 


人 
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波束 斌 形 
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pi >=0 ;|u|=1; pr 20;pry <1;p,>0 


SINR,,, 是 PU1 的 虚拟 SINR。 由 pi, 的 特性 可 知 , 问题 P2 
和 P3 有 相同 的 结果 。 因 为 当 p=1 时 ，P3 即 为 P2。 为 确保 
px =1， 在 算法 中 增加 |p,,, 一 了 | 的 收敛 性 判 据 。 
1.3 基于 上 行 和 下 行 对 偶 的 问题 转换 

根据 上 行 和 下 行 的 对 偶 特性 60, 依据 文献 [7] 采 用 上 行 虚 拟 
功率 ， 下 行 问题 P3 可 转换 为 虚拟 上 行 问题 P4。 


. K ， L 
P4: 六 7poPo + Oi ) qi - DPruges (9) 
qu Ru 
Ss.t. SINR! = K+ >1; k=1,2,...,K 
| > qrB, + 
1 天 
gq: =>0;|u|=1; pr, >0; pr <1; qr >0;p, >0;1=1,2,...,L 


gag 和 qu sgk 是 分 别 对 应 于 个 SU 和 工 个 PU 的 
虚拟 上 行 功率 。 

不 同 于 文献 [7] 的 是 ， 
的 nn 变 为 


于 考虑 了 PN 对 SU 的 影响 ， 其 中 


1 三 [CPP to) NPR + oD rk=1,2,..,K。 


了 关 优 化 问题 P4 的 两 点 说 明 : 
a)P4 和 P3 有 相同 的 优化 波束 赋 形 向 量 w (k=1,2,...,K )， 
通过 P4 可 以 较 简便 地 得 到 。 
b)P3 的 P (CCK=12,R ) 和 pp， 《1=1,2,..,L) 由 如 下 的 
关系 式 得 到 : 
设 p=(p,py.wpx) 是 个 SU 的 下 行 功 率 向 量 ; 
pp =(pispisw pes) 是 上 个 PU 的 虚拟 下 行 功率 向 量 , 有 


p=(D-G)'n (10) 
b=G,p, (1D 
其 中 : [DY]. 一 ur Ru ; [G,]; yeuU’ RU (12) 
. i=j 
= 13 
[G]， bs Ru i#j (13) 
n=[yCp,R to),.s Np,R + od) (14) 


[=1,2,..,L ,ki,j=1,2,...,K, 
2 ”和 迭代 算法 
2.1 算法 思路 
P4 可 知 ， 当 q, 固定 时 ， 波 束 赋 形 向 量 以 的 值 应 当 使 
SIVRY 尽 可 能 的 大 ; 而 当 u, 固定 时 ，g, 将 在 尽量 小 的 情形 下 满 
足 SINR" >1 的 要 求 。 SINR! 的 最 低 要 求 是 等 于 1， 因此 可 得 到 
由 和 和 《k=1,2,…, ) 的 迭代 算式 。 而 gr 是 因 六 ,所 生成 ， 
所 以 只 需 mvgx， 即 可 生成 迭代 算式 。 
2.2 ”算法 流程 
迭代 算法 如 下 : 
a) 初始 化 gq (0)=1，n=12,...,K+L 
b) 迭代 次 数 z=12,…， 作 以 下 的 和 迭代 步 又 ， 直 到 按 要 求 收 


lm 
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(a) 波束 赋 


NY 
Ba 


子 更 新 : 


qi(t—Dur DR CD 


一 InaX +L a 
Hi(D ux D1 “ol py q(t-DyR 四 ;kl;2,..K 
i=1,iz#k 
(15) 
对 于 上 式 求解 ， 由 文献 [10, 14] 可 知 。 4.(D) 和 ui0) 是 矩阵 


对 区 二 了 IC-D7ZP 圳 的 最 大 广义 特征 值 和 对 应 


的 单位 特征 向 量 。 
(b)》 SU 的 虚拟 上 行 功率 分 配 更 新 : 


on (16) 
(c) 转换 为 下 行 功率 分 配 : 
PD)=(DD -Gn (17 
Pb,(D)=G,(D PD (18) 
其 中 : 
[DO = Ue WRUD ;GD = yxy Rx,u, CD) ;(19) 
aon, = er. 本 C0) 
1=1,2,.., Lk,i,j=1,2,..,K 
(d) PU 的 虚拟 上 行 功率 分 配 更 新 : 
A 
c) 收敛 性 检测 。 
[qi -g(t -Dz<6;|pri) -1<6 (22) 
=1,2,.., 民 + 工 ;1=1,2,..…, 工 ;6 是 一 预先 定义 的 收 钱 判决 门限 。 


若 满足 收敛 条 件 ,到 ;否则 , 转 到 b)。 
d) 当 经 过 N 次 透 代 后 , 满足 收敛 条 件 时 ， 可 得 问题 最 优 值 


为 : 


(a) 最 优 功率 分 配 : jp (N),p,(N),…, Px(N); 
(b》 对 应 的 波束 赋 形 向 量 vy (N),w(N),.…,Ux(N)。 
2.3 ”算法 收敛 特性 分 析 


定理 ”算法 收敛 需要 满足 条 件 lim j4 (7) 一 1; 


k=],2v65 RK 


和 lim py ->1;， 71=1 工 。 


证 明 由 了 (D= 并 (全 i ,sR ) 可 知 : 
Mi (DD 
lim MD >1 = lim a | 
1 一 oo 1 


JL(D) 是 实时 的 SINR? ,和 途 代 过 程 中 的 收银 过 程 如 下 : 
ad-DU (Ru CD 

= Max K+L 可 知 
HD Tao “ol py -dR tu wy. 


i=1,i#k 


A 是 q, 的 单调 递增 函数 ， 是 g (zk,i=1,2,.…,K+D) 的 单 
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季 中 性 
调 递减 函数 。 由 及 的 性 质 ， 可 知 J(n) >0。 在 迭代 的 第 ; 步 ， 
JL) 的 取 值 有 两 种 情形 ， 分 别 是 
a) AD >1 
由 了 (= 和 DD 12,.,K， 可 知 : 


oA) 


gS q(t-D k=12,...,K 


对 于 gi ,0); 1=1..,L， 由 迭代 关系 : 
二 Pru(Dqru(t—b) if min{p,..., pr} >0 
dk+1 (0D) 。 
dk 一 了 otherwise 


可 知 : gx, (0) < gry(t-D 或 gx,() > qrult-D 
第 1+1 步 时 ， jy.(t+1) 有 如 下 的 两 种 情形 : (a) jy (t+1) < 
如 0， 可 得 ID 减 小 ，(b) (t+DD)>p(1)， 可 得 gq() 仍 是 减 


小 。 

情形 (a) 的 j 将 由 于 gq, 的 减 小 趋 于 1 (b) 的 j 将 由 于 
qi 的 减 小 在 迭代 若干 次 后 改变 不 等 号 而 成 为 情形 (a)， 否 则 ， 
当 g 一 0 时 ， 有 人 一 0， 与 a) 的 前 提 条 件 矛 盾 。 

b) 0<y0) <1 

可 和 gD)> gt-D); k=1,2,...,K 

第 1+1 步 时 ，w+D 有 如 下 的 两 种 情形 : 


a) 人 (+]) > 


L(t) ， 可 得 gq (7) 是 增加 的 ; b) w+D<(D ， 可 得 g(D 仍 
是 增加 的 。 
在 迭代 过 程 中 , 情形 (a) 的 jy 将 由 gq 的 增加 趋 于 1; (b) 


的 六 将 由 于 4 的 增加 在 迭代 若干 次 后 改变 不 等 号 而 成 为 (a)， 
否则 ， 当 gw 时 ， 有 j 与 b) 的 前 提 条 件 矛 盾 。 
二 人 if min{ pe}>0 
qru(t—D) 
当 p>1 时 ， 表 明 PU 所 受到 的 干扰 增 大 ， 需 要 减 小 p, ;jp 
的 减 小 致使 p, ,也 减 小 。 当 0< p<1 时 ， 表 明 PU 所 受到 的 
干扰 减 小 , 应 当 增 加 p, 来 得 到 较 高 的 SINR, ; 由 此 , 相应 的 六， 
也 增加 。 


> oo ， 


可 知 ， 


otherwise 


I 


lim Pru(D 一 1 时 ， lim [qxu(t-D-gqrn D0 ; 


1=1,2,..,L 
所 以 ， 算 法 收敛 。 
证 明 完 成 。 
2.4 算法 复杂 度 分 析 
以 乘 和 加 作为 基本 运算 单位 ， 算 法 单 次 迭代 过 程 的 各 步 又 
的 计算 复杂 度 03 如 下 : 
正 值 和 


a) 区 q(t—DYyR 圳 的 最 大 广义 特 和 


对 应 的 特征 向 量 求解 。 ne 其 中 运算 量 
最 大 的 是 N, x N 矩阵 求 道 ， 因 此 运算 复杂 度 是 O(N3) 。 当 


k=1,2,.…K 时 ， 运 算 复杂 度 是 O(KN;)。 
t—1 次 运 
) = 全 和 -1,2,..,k， 共 作 色 次 运算 。 
Hi (D) 


录用 稿 


c) [D4=Us (DRu(),， 共 作 2 KN?1+2 KN, -次 运算 。 
d) [GODR LUGOD :1 共 作 2 LKN? +2 LKN, 次 


运算 。 


e) [GOD),, he 0 
NU (DRu,D, 
-2 KN? -2 KN, 次 运算 。 
f) p(D)=(D(WD -GD) in， 共 作 3+2 kK? 次 运算 。 
g) 记 (D) =G,(Dp(r)， 共 作 2 LK 次 运算 。 
h) 0 if min{ ,Prx}>0 
qru(t—Db) 


共 作 2 K?N? +2 KN, 


， 共 做 


otherwise 
L+ 久 次 运算 。 

i) 收敛 性 检测 ， 共 做 2 g +4 上 次 运算 。 

因此 ， 在 一 个 迭代 过 程 中 ， 算 法 步骤 b)~5D 的 运算 次 数 是 2 
天 2N 12 LKN? +2 K°N, +2 LKN, + 


Kt2K*+2LK1SLT3K. 
通常 情况 下 ， 久 和 的 数量 级 相同 ， ~、>1， 可 以 略 去 
K?N,、LKN,、K3、K?、LK、 工 和 K 等 低 阶 项 , 并 
可 得 算法 步骤 b)~D 的 复杂 度 是 O(K?N?)。 


了 二 天， 


季 中 恒 ， 等 


， 基于 上 下 行 对 偶 的 CR 下 行 劝 率 
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ZX 


a WN 


法 有 可 行 解 的 信道 实现 所 占 百分比 。 它 是 算法 对 信道 变化 适应 


性 的 一 种 度量 。 


取 y,=y 


Feasibility Percentage 
o 
a 


(=1,2,...,K )。 图 中 的 gamma 即 为 


ee 
—e— Pp- 
nie 


30dBm 
27dBm | 
0 


图 2 yy 与 可 解 度 百分比 的 关系 
如 图 2 所 示 ， 当 PBS 的 发 送 功率 变化 时 ， 可 行 解 区 域 也 发 
生变 化 。 这 是 由 于 PBS 的 发 送 功率 p, 的 引入 ， 影 响 了 优化 问 
题 的 信 干 噪 比 , 进一步 影响 了 CRN 的 干扰 约束 , 使 可 行 解 区 域 
缩小 .如 以 图 中 的 可 解 度 百分比 为 0.6 为 例 , 当 p, =30 dBm 时 ， 
相应 的 y 约 为 4.0; 而 当 


步骤 a) 的 复杂 度 O(KN;) 可 得 算法 的 复杂 度 是 
O(max(KN3,K?N?)) 。 而 由 文献 [7] 知 ，SDR 算法 的 复杂 度 是 
OCK3N3) 。 因 此 ， 和 迭代 算法 的 复杂 度 要 比 SDR 算法 低 很 多 。 


3 ”数值 仿真 及 结果 分 析 


数值 仿真 条 件 : SU 用 户 数 及 =5， 由 =5; CBS 的 


发 送 天 线 数 N, =9; 


PU | 


户 的 干扰 门限 =-6 dB ， 


K+l 


(C1=12,, 工 ); SU 的 QoS 约束 y, (k=1,2,.….,K ) 在 区 间 [0.1， 


9.0] 内 取 值 ， 噪 声 功 率 o? =-6 dB，( k=1,2,.…,K); 信道 是 
Rayleign 衰落 信道 ， 对 于 每 一 次 的 信道 计算 结果 ， 均 作 100 次 
独立 Monte-Carlo 运算 再 取 均 值 .PBS 的 发 送 功率 p, =30 dBm。 
算法 的 收敛 门限 5=104。 以 下 图 3 中 的 The sum of P1 即 是 
标 值 六 。 
全 

3.1 和 迭代 算法 的 收敛 性 

迭 代 算 法 的 收敛 门限 和 友 代 次 数 如 表 1 所 示 。 

表 1 不 同 收敛 门限 下 的 迭代 次 数 
类 别 组 1 组 2 组 3 组 4 组 5 组 6 组 7 
收敛 门限 8 102 103 104 105 106 107 10 
迭代 次 数 N 9 14 19 34 36 58 6 


表 1 中 的 各 组 数据 可 知 ， 算 法 收敛 很 快 。 如 当 收 敛 门限 
为 104 时 ,经 过 36 次 迭代 即 满足 收敛 条 件 。, 相 对 于 SDR 算法 ， 
迭代 算法 不 仅 降低 了 计算 量 ， 而 且 收 敛 速 度 快 。 
3.2 PBS 的 发 送 功率 与 可 行 解 区 域 的 关系 
可 行 解 区 域 的 度量 指标 是 可 解 度 百 分 比 〈feasibility 
percentage)， 是 指 在 一 定数 量 的 随机 信道 实现 条 件 下 ， 优 化 算 


ul 


在 无 PBS 时 (p,=0), y 


( p, =30 dBm 和 p, =27 dBm)， 可 解 度 百 分 
要 获得 y =7.0 的 信 干 品 比 ， 


大 导致 CBS 对 PU 的 干扰 


泥 可 
解 区 域 的 影响 是 
3.3 ”PBS 的 发 送 

如 图 3 所 示 ， 当 
功率 增加 。 


PBS 对 SU 的 干扰 也 逐 
增加 下 行 功率 。 这 表明 在 
的 QoS 有 明显 的 影响 。 


Ps =27 dBm 时 ,相应 的 y 增加 到 约 4.7; 
能 到 约 7.0。 而 对 于 存在 PBS 的 情形 
比 为 0。 这 是 因为 
功率 必须 增加 ， 而 下 行 功率 增 
增 大 ， 无 法 满足 干扰 约束 条 件 ， 所 以 
行 解 ， 可 解 度 百 分 比 为 0。 由 此 可 见 ，PBS 的 功率 对 可 行 
很 显著 的 。 
功率 与 下 行 
PBS 的 发 送 功率 p, 增加 时 ，CRN 的 下 行 
当 p, <30dBm 时 ， 下 行 功率 增 力 
>30dBm 后 ， 下 行 功率 增加 较 快 。 这 是 因为 当 
簿 渐 增 加 ，CRN 要 保证 


下 行 


功率 的 关系 


上 较 平 缓 ; 当 


p, 
p, 逐渐 增 大 时 ， 


SU 的 SINR， 必 须 


CRN 系统 中 ，PBS 的 


发 送 功率 对 SU 


The sum of P1 


图 3 


4 ”结束 语 


过 对 工 


30 
Pp(dBm) 


Ps 与 下 行 功率 的 关系 


赋 形 问题 的 研究 ， 得 到 了 


PN 干扰 条 件 


便 、 快速 和 实 


工具 ， 利 于 工程 实现 。 针 对 PN 的 干 # 


j 的 迭代 算 泣 


作 于 underlay 模式 的 CRN 下 行 功率 分 配 和 波束 
的 求解 优化 问题 
,为 CRN 的 组 网 提供 了 一 种 有 效 的 


的 简 


效应 ， 


进行 了 算法 数值 


录用 稿 


仿真 , 得 到 了 PBS 的 发 送 功率 变化 对 可 行 解 区 域 和 下 行 功率 均 
有 较 明显 的 影响 的 实用 结论 ,进一步 表明 PN 的 干扰 在 CRN 网 
络 模型 中 是 不 能 忽略 的 , 尤其 是 在 CRN 的 网 络 构建 工程 中 , 更 
应 当 重 视 。 
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